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The analysis of the 21 cm signature of cosmic string wakes is extended in several ways. First we 
consider the constraints on Gfj, from the absorption signal of shock heated wakes laid down much 
later than matter radiation equality. Secondly we analyze the signal of diffuse wake, that is those 
wakes in which there is a baryon overdensity but which have not shock heated. Finally we compare 
the size of these signals to the expected thermal noise per pixel which dominates over the background 
cosmic gas brightness temperature and find that the cosmic string signal will exceed the thermal 
noise of an individual pixel in the Square Kilometre Array for string tensions Gfi > 2.5 x 10~ 8 . 
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PACS numbers: 98.80.Cq 



I. INTRODUCTION 



^C) In a previous paper pQ we studied the signal of a shock heated cosmic string wake in a position space 21cm redshift 

maps. We pointed out that cosmic string wakes give rise to a pronounced signal in such maps, in particular at redshifts 
larger than that of conventional reionization. The signal of a single wake is a wedge which is extended (a degree or 
larger for redshifts of z ~ 20) in both angular directions and thin in redshift direction. This effect comes from the 
fact that wakes form overdensities in baryons already at very high redshifts, and that these regions then lead to the 
extra 21cm emission or absorption. 

In a followup paper [5] we computed the angular power spectrum for shock heated cosmic string wakes emitting 21 
cm radiation. Because of our interest in shock heating in both [T] and [3] we focused our attention on those wakes with 
the greatest kinetic gas temperature Tk- Thus we concentrated our analyses on the wakes formed earliest, Zi ~ 3000 
and emitting 21 cm radiation as late as possible before reionization, z e ~ 20. This led us to miss the constraints on 
G[i that shock heated wakes formed at later Zi could provide by considering the absorption of 21 cm radiation. In 
— ' this paper we consider these constraints. 

Furthermore, even if baryons collapsing onto the primordial string wake do not shock heat, the string will neverthe- 
less induce an overdense region of baryons. This region, however, will be more diffuse than if shock heating occurs. 
The wakes will be thicker but the overdensity smaller. In this paper, we study the 21cm signal of such a "diffuse" 
C^ cosmic string wake. 

Current constraints on the cosmic string tension from analyses of the angular power spectrum of CMB anisotropy 
CN maps yield 

Q Gfi < 1.5 x 10~ 7 (1) 



using combined data from WMAP and SPT microwave experiments [3J. Limits based only on WMAP data gives a 
bound larger than this by a factor of two [4] . It must be emphasized that any limits on G/i coming from power spectra 
analyses implicitly depend on parameters describing the cosmic string scaling solution which are quite uncertain. 
Two prime examples of such parameters are the number N of string segments crossing any Hubble volume, and the 
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curvature radius c\ of a long string relative to the Hubble radius. These parameters are known only to within an 
order of magnitude since they must be determined via numerical simulations of cosmic string network evolution, and 
these simulations are extremely challenging because of the huge hierarchy of scales which must be included at the 
same time. 

In our analyses below we show that a cosmic string signal will exceed the thermal noise of an individual pixel in a 
future radio telescope such as the Square Kilometre Array (SKA) for a string tension 

G^>2.5xl0" 8 . (2) 

Our new limit is independent of the parameters N and c\. In addition, since the signal of a string wake has a very 
special geometry in position space, position space shape detection algorithms are likely to be able to detect string 
wake signals even if they do not stick out above the Gaussian noise on a pixel by pixel basis. The situation is similar 
to what is encountered when searching for cosmic string signals in CMB anisotropy maps: since long string segments 
produce line discontinuities in the maps, these can be identified by edge detection algorithms such as the Canny 
algorithm [5 for values of the tension substantially lower than the value for which the string signal dominates on a 
pixel by pixel basis, as has been studied in detail in [6 a . 



The paper is organized as follows: We begin in Section ITT] with a brief review of cosmic strings, their wakes, and 
their 21 cm signatures. Readers familiar with such results may skip this section. Section |III| is the main part of our 
work. In this section we extend the analysis of the 21 cm signature of cosmic string wakes in several ways. First we 
consider the constraints on G/j, from the absorptions signal of shock heated wakes laid down much later than matter 
radiation equality. Secondly we analyze the signal of diffuse wakes, that is those wakes in which there is a baryon 
overdensity but which have not shock heated. Finally we compare the size of these signals to the expected thermal 
noise per pixel which dominates over the background cosmic gas brightness temperature. In Section |IV| we present 
our conclusions and put our work in context. 

II. COSMIC STRINGS, THEIR WAKES AND 21 CM SIGNATURES 

Many particle physics models beyond the Standard Model predict the existence of cosmic strings. In particular, 
cosmic strings are produced at the end of inflation in many models of branc inflation [7] and in many supergravity 
models [5]. Cosmic strings may also survive in early universe models based on superstring theory such as "string 
gas cosmology" UH] • Since the amplitude of cosmological signatures of cosmic strings is proportional to the string 
tension, which itself is set by the energy scale of the new physics, searching for signatures of strings in cosmology is a 
way of probing Beyond the Standard Model physics which is complementary to accelerator probes such as the Large 
Hadron Collider experiments which probe the low energy limit of such theories. 

Cosmic strings can be used as a cosmological probe of New Physics because any particle physics model which 
admits stable cosmic strings as a solution inevitably results in a a network of such strings forming during a symmetry 
breaking phase transition in the early universe (see |11H13| for reviews on cosmic strings and cosmology). Such a 
network will persist to the present time and approaches a dynamical attractor configuration, the so-called "scaling 
solution" in which the statistical properties of the string network are independent of time if all lengths are scaled to 
the Hubble radius. The cosmic string scaling solution is charaterized by a random-walk-like network of infinite (or 
"long") strings with step length comparable to the Hubble radius, plus a distribution of cosmic string loops with radii 
smaller than the Hubble radius. The typical transverse velocity of a long string segment is the speed of light. This 
implies that at all times t > t eq (t eq being the time of equal matter and radiation) relevant for structure formation 
there are strings which act as seeds for the accretion of matter. The key parameter which characterizes cosmic strings 
and their cosmological consequences is the mass per unit length fi. It is usually quoted in terms of the dimensionless 
number Gfi, where G is Newton's gravitational constant. 

Both long string segments and string loops lead to nonlinearities in the matter distribution at all redshifts up to z eq , 
the redshift of equal matter and radiation. Cosmic string network simulations 14J indicate that the long strings have a 
larger effect than the loops, so we will here focus our attention on the effects of long string segments, as we did in [TJH] 
(For a recent study of the 21cm signal of a cosmic string loop see [IS]-)- -"- n the standard ACDM cosmology without 
strings, nonlinear structures are negligible until shortly before the time of reionization. Hence, the signals of cosmic 
strings should be more and more pronounced against the "noise" of effects from structure formation in the ACDM 
model the higher the redshift is. 21cm surveys hence appear as an ideal window to probe for the possible existence of 
cosmic strings. They probe the distribution of matter at much higher redshifts than galaxy redshift surveys since they 
map out the distribution of neutral hydrogen in the universe, in particular before the time of reionization. Compared 
to CMB maps, 21cm maps have the advantage of yielding three-dimensional information. 

The nonlinearity in the matter distribution induced by a long string segment at time t takes the form of a "wake" 
[16] whose length is set by the length of the string segment (which will be similar to the Hubble radius t), whose 
width is set by v s "/ s t, where v s is the transverse velocity of the string segment and j s is the associated relativistic 
gamma factor, and whose average thickness is proportional to Gfit. The wake formation starts with an initial velocity 
perturbation and then grows by gravitational accretion. 

The velocity perturbation (see Figure [l]) stems from the fact that space perpendicular to a long string segment 
is conical, i.e. flat space with a missing wedge whose deficit angle is a = 8ttG/j,. A string moving through a gas of 
particles with transverse velocity v s will hence induce a velocity perturbation Sv of the gas towards the plane behind 
the string, whose magnitude is given by Sv = 4irG[xv s j s . . 

A cosmic string segment laid down at time ti (we are interested in ti > t eq ) will generate a wake with physical 
dimensions: 

h(U) x l 2 (ti) x w(U) = U ci x ti v s ~f s x U A7rGiJ,v s ^ s . (3) 

where c\ is a constant of order one. In the above, the first dimension is the length in direction of the string, the 
second is the depth, and the third is the average width. After being laid down, the wake will grow by gravitational 
accretion. This has been studied using the Zel'dovich approximation |17j in a number of works |18j . In this method 
of analysis, we follow the time evolution of the height above the centre of the wake of mass shells to determine the 





FIG. 1: Sketch of the mechanism by which a wake behind a moving string is generated. Consider a string perpendicular to 
the plane of the graph moving straight downward. From the point of view of the frame in which the string is at rest, matter 
is moving upwards, as indicated with the arrows in the left panel. From the point of view of an observer sitting behind the 
string (relative to the string motion) matter flowing past the string receives a velocity kick towards the plane determined by 
the direction of the string and the velocity vector (right panel). This velocity kick towards the plane leads to a wedge-shaped 
region behind the string with twice the background density (the shaded region in the right panel). 



comoving distance q n i which "turns around" at time t, i.e. for which the physical height h(q, t) is maximal. We obtain 
that the comoving thickness grows linearly with the scale factor of the universe. A wake laid down at redshift Zi will 
have comoving thickness at redshift z given by 



16tt Gfiv s j s 
q n i{z) = — — (Zi + l) 
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(4) 



The physical height of the shell above the centre of the wake at the time of "turnaround" is half what the height would 
be in the absence of gravitational accretion, and hence the relative overdensity is a factor of 2. It is then assumed 
that the shell collapses to half of the height it had at turnaround, that it then undergoes shocks and virializes at 
that distance. Then, the overdensity will be a factor of 4 greater than the background density. This approximate 
analytical analysis has been confirmed in detailed hydrodynamical simulations |19j . 

During the infall onto the string wake, the baryons acquire kinetic energy which during the process of shock heating 
transforms to thermal energy. The kinetic temperature T% acquired by the hydrogen atoms can be obtained from 
the kinetic energy which particles acquire when they hit the shock front. This is determined using the Zel'dovich 
approximation IT]. The result obtained in pQ is 



T K ~ [20K](G M )2( Us7s 



,Zj + 1 

z + 1 



(5) 



where (Gfj,)e indicates the value of G/i in units of 10 -6 . The result is expressed in degrees Kelvin. 

Note that the wake kinetic temperature increases as z decreases since the wakes have more time to grow, which 
in turn leads to a growing gravitational potential. Since T 7 decreases in time, then it is more and more likely that 
the string wake 21cm signal will be in emission rather than absorption as t increases. However, it turns out that 
for reshifts significantly larger than that of conventional reionization and for string tensions smaller than the current 
limit (TH), Tk is smaller than T 1 and that hence the string wake signal will be in absorption. 

The 21cm signal of a cosmic string wake has a distinctive geometry. We will see the string-induced absorption (or 
emission) regions in directions of the sky where our past light cone intersects the cosmic string wake. These are wedges 
in the sky of angular scale determined by the comoving Hubble radius at the time that the string is formed. The 
string wake grows in thickness, but its planar dimensions are constant in comoving coordinates. For wakes formed at 
a redshift corresponding to recombination this angular scale is about one degree, for wakes formed at equal matter 
and radiation slightly smaller. Since photons emitted at the top and bottom of the wake undergo slightly different 
redshifts, the region in 21cm redshift surveys of extra emission/absorption takes on the form of a thin wedge in redshift 
direction with a thickness which is [1] 
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III. DIFFUSE COSMIC STRING WAKES 

The analysis in the previous section neglected the temperature of the gas at the time that the wake was formed. 
The accretion onto a string wake further increases the intrinsic gas velocities, leading to an effective gas temperature 
which is 2.5T gi where T g is the temperature of the background gas. If 

T K < 2.5 x T g , (7) 

then the incoherent velocities due to the thermal motion of the accreted gas dominate over the coherent velocity 
induced by accretion. In this case, there is no shock heating, and the thickness of the overdense region induced by 
the wake is larger than it would be in the absence of thermal motion of the gas particles. The resulting overdense 
region we will call a "diffuse wake" . 

We now estimate the size of the diffuse wake. In linear cosmological perturbation theory, the total mass accreted 
by the string wake is the same both in the cases with and without shock heating. The difference is that the large 
intrinsic thermal gas velocity will render the wake wider and hence less dense. 

The first step is to compute the width h w of the wake obtained taking into account the gas temperature T g . The 
result is 

h w (z)\T K <T g = h w {z)\ Tg=Q (|l) (8) 

for Tk < T g . The diffuse wake will thus be thicker by the ratio of temperatures T g /Tx, but the density of baryons 
will be smaller by the same factor. Hence, the width of the 21cm signal induced by a diffuse wake will be thicker in 
redshift direction by the above ratio of temperatures compared to what is given in ([6]). On the other hand, the 21cm 
brightness temperature will be smaller since it is proportional to the baryon density which is smaller than in the case 
of a shock-heated wake. 

The formula Q for the width of a diffuse wake can be derived by assuming equipartition of energy between thermal 
energy and potential energy which for Tk <C T g yields 

m H i8* = \t 9 , (9) 

where mui is the mass of a hydrogen atom and <5<£> is the gravitational potential induced by the wake. The latter, in 
turn, is 

(5$ = 2irGa\h\ , (10) 



where a is the surface density of the wake, and \h\ is the height. Combining (pjl) with (10) yields the linear scaling in 
temperature given in (J8J). 

The surface density a of a diffuse wake is the same as that of a shock heated wake. It is given by the energy per unit 
area of matter within a comoving height q n i(z), the shell which is turning around at redshift z. Thus <j{z) — q n i{z)po, 
where p is the current background density. The overdensity Ap in the diffuse wake is given by 



which results in 



'-'l + £'). (12) 
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The extra baryon density in the wake leads to extra 21 cm absorption or emission. 

The expression for the wake brightness temperature written in terms of kinetic temperature Tk s and the hydrogen 
atom number density tihi is very similar for a shock heated or a diffuse wake. The only difference is that Tx g is 
the kinetic temperature Tk of the wake when shock heating occurs whereas it should be interpreted as the kinetic 
temperature of the cosmic gas T g for diffuse wakes. 

T 7 \ n H f e (1 + Ze) 1 / 2 



ST b (z e ) = [17 mK] — 1 2_ -^- y e > — . (13) 

V ' [ i l + x c \ T Kg ) n b B T 2sin 2 r 
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Here, z e is the redshift of 21 cm emission or absorption, T 7 is the CMB temperature, 9 is the angle of the 21 cm 
ray with respect to the vertical to the wake, and x c are the collision coefficients whose values can be obtained from 



the tables listed in [20]. Here and throughout we take the cosmological parameters to be Hq = 73 km s~ Mpc -1 , 
ilb = 0.0425, O m = 0.26. We work with a matter dominated universe for z < 3000 with the age of the universe 
to = 4.3 x 10 17 s. The origin of the sin~ 2 (#) factor will be discussed in the Appendix, where it will also be shown that 
this factor does not lead to any physical divergence. 

For a shock heated wake, the ratio of the wake's hydrogen atom number density to the background density, UHi/n^r, 
is 4. For a diffuse wake the relative brightness temperature is lower than for a wake which has undergone shock heating 



background density njij is smaller. For our analysis below we take: 



since the number density tihi of hydrogen atoms inside the diffuse wake is smaller, and thus its ratio in ( 13 ) to the 
bg 
l HI 



n W=(^^) T K <3T g (14) 

n HI 9 

= 4 T K > 3T g (15) 

Note that for simplicity we have used T% < 3T„ instead of Tr- < 2.5T„. 



The brightness temperature signal in ( 13 ) discussed in [TJ [2] was compared to the average brightness temperature 
of the surrounding cosmic gas. Between redshift z — 20 and 35 the average background brightness temperature varies 
from -0.34 mK to -8.6 mK, respectively. However the analysis of the signal did not consider beam size nor thermal 
noise. In fact for z ~ 20 it is the thermal noise, and not the background brightness temperature, that will first limit 
our ability to detect a cosmic string wake signal. The opposite is true for z = 35. Beam size and thermal noise for 
the 21 cm brightness temperature were considered in |21j where their Eq. (4.19) gives the thermal noise per pixel: 



T„ 



12 mK (arcmin/6> rcso i u tion) fl + z 



y/(T/Whr){A e /km 2 ) V 21 J ly/l + z-1 



3.85 r 



VT 



1/2 



(16) 



Here r is the total observing time and A e is the effective antenna area. To derive this noise temperature per pixel 
in [21] we began with the thermal noise per visibility as given by Morales [22] to arrive at the thermal noise per 
pixel for the brightness temperature. The angular resolution is X/D where A is the observation frequency and D is 
the baseline. The angular resolution is assumed to be tuned by a dilution of the array from being fully compact by 
a simple scaling of all baseline positions by a fixed amount. The system temperature is given by ARCADE 2 |23j . 
Further details can be found in [21] . 

The noise per pixel resolution dependence given in eq. [16] assumes an isotropic resolution in all three spatial 
directions. However for a cosmic wake signal, the width is order AnG^Zi + l)/(z + 1) smaller than the two length 
directions, whereas each length direction is of order the Hubble size. If 9hubb\c(z) is the angular resolution needed 
for the wake's length, then #hubbie(z)47rG /i(zj + l)/(z + 1) is the angular resolution needed for the width. Thus an 
appropriate estimate of the resolution in radians needed to calculate the noise per pixel is: 

/ 47rG>(z i + lh i/3 / 47rg^+l) y/3 1 

^resolutions - { {z + ^ ) *hubbto(*) - [ £ + 1} J (,/F + T - 1) ( 17 ) 

Hence, in a three dimensional map with the necessary resolution to detect a cosmic string wake laid down at Zi and 
emitting or absorbing radiation at z, the thermal noise per pixel is: 

*- - *^Wm^^r* - w^ 

where (G^) 9 = 10 9 G>. 

We will consider 10 000 hours of total observing time. Evidently for the SKA the effective antenna area is 1 sq 
km. For Zi — 3000 and z — 20, a G[i = 4 x 10~ 8 leads to an average noise per pixel T n — 0.3 mK comparable in 
size to the background brightness temperature. Smaller G/i require a smaller resolution and lead to larger noise per 
pixel. Hence for G/j, < 4 x 10~ 8 the thermal noise dominates over the background brightness temperature for the 
wake signal. For example at Zi = 3000 and z — 20 a G[i = 5x 10~ 9 gives a ^resolution (z) — 5.8 x 10~ 3 radians, which 
is 20 minutes of arc, and hence an average noise per pixel of T„ = 0.66 mK. For strings laid down at a lower redshift 
the noise is greater. The noise is also a rising function of the emission redshift z. For a Zi — 1000 with z — 20, 25 and 
35, T n = 1.6 mK, 2.0 mK and 2.7 mK, respectively. 

In Figured we show the relative brightness temperature (vertical axis) as a function of (G/j,)q (horizontal axis) for 
various values of the redshift Zi of wake formation, evaluated at the redshift z — 20. Negative (positive) brightness 
temperature means absorption (emission). The two almost horizontal (brown) lines give the noise per pixel level in 
an experiment such as the SKA for the pixel size chosen such as to optimize the search for a string wake with the 
respective value of Gfi (see the above discussion). 
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FIG. 2: The relative brightness temperature (vertical axis) in degrees Kelvin as a function of (Gfi)e, (horizontal axis) for various 
values of the formation redshift Zi, evaluated for an observation redshift of z = 20. The curves from left to right (in the region 
of low values of G/x) correspond to z% — 3000 (red curve), Zi = 2500 (orange), Zi = 2000 (green), Zi = 1500 (light blue) and 
Zi — 1000 (dark blue). The two brown lines indicate the expected thermal noise per pixel in an experiment such as the SKA, 
with a pixel size chosen to depend on G[i as described in the text. The black line that is almost indistinguishable from the 
x-axis is the brightness temperature of the background cosmic gas. 
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FIG. 3: A zoom in of figure [2] to the point where the noise curve (brown) crosses the signal curves (colours red to blue). The 
black line at the top of the graph is the brightness temperature of the background cosmic gas. 



As expected, the larger the value of Zi, the easier it is to detect low values of Gfi since the wakes have had a 
longer time to grow in thickness. A feature overlooked in our previous work [TJ [2] is that for a fixed observational 
redshift z, there will be a range of values of G[i for which the wake's brightness temperature Tjc will be so close to the 
background cosmic gasp temperature that there will be no signal. In this situation, we can consider smaller values of 
Zi which lead to a clean and large absorption signal. 

For each curve in Figure ^ it is apparent that there is a kink in the brightness curve (in the absorption region) . 
This kink occurs at the value of G/x for which Tk — 3T ff , the transition point between a wake with shock heating 
(larger values of Gfi) and a diffuse wake. Our previous analysis was only valid until that point. We see here that with 
our current work it is possible to extend the range of values of Gfi which can be probed by 21cm redshift surveys. 

Taking into account that the optimal pixel size will depend on Gfi, then the pixel noise level will also depend on 
Gfi. In Figure [3] we zoom in on the small Gfi region of the previous figure and compare the string signal to the pixel 
noise. We see that the absorption signal of the diffuse wakes is larger than the noise level for values of Gfi larger than 



Gfi > 2.5 x 10" 



(19) 



for a detection redshift of z — 20. This is an improvement of the current bound M and demonstrates the potential 
of 21cm surveys to probe cosmic strings. 



For smaller values of Gfj,, cosmic string wakes will still provide a signal. The signal has a particular shape in a 
three-dimensional redshift survey and should hence be detectable using shape detection algorithms for values of Gfi 
much smaller than the above value (19). In work in progress we are exploring this issue quantitatively. Note that a 



similar problem arises in the case of identifying cosmic string signatures in CMB temperature maps. It was shown 
[6] that edge detection algorithms provide the promise of pushing the limits of detection several orders of magnitude 
below the level of the pixel noise, and one order of magnitude lower than limits which can be achieved using CMB 
angular power spectra studies. Thus, we are optimistic that 21cm redshift surveys will provide a very promising arena 
to probe for the possible existence of cosmic strings. 
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FIG. 4: The relative brightness temperature (vertical axis) as a function of formation redshift Zi (horizontal axis) for various 
values of (G/J,)e,- The values of (Gfi}6 chosen are 0.05, 0.09, 0.13, 0.17, 0.21, and 0.25 (red, orange, green, light blue, dark 
blue and purple curves, respectively). At low values of z% some of these curves are labelled by an increasing amplitude of the 
absorption signal. The black and brown lines near the x-axis represents the background gas brightness temperature and the 
noise temperature, respectively. They are nearly indistinguishable on this scale for these values of G[s,. The noise is plotted for 
(Gfi)e = 0.05. The larger values of (Gfi)e would give less noise. 

In Figure |4] we show how that formation redshift which leads to the largest amplitude of the absorption signal 
depends on Zj. If we are looking for strings of a particular value of G/j, we would be focusing on wakes formed at 
the redshift which yields largest absorption amplitude. Since the formation redshift determines the size of the 21cm 
wedge which the corresponding wake produces, we can use the size information to optimize the search strategy. 



IV. CONCLUSIONS 



We have extended the analysis of the 21cm signal of cosmic string wakes to the case when the string tension is too 
small to yield shock heating of the wake. Instead, a "diffuse wake" forms. This wake is less dense but thicker than a 
wake which has undergone shock heating. Hence, the relative 21cm brightness temperature is smaller in amplitude, 
but it is extended over a larger redshift interval. 

Demanding that the relative brightness temperature exceeds the pixel thermal noise in a future 21cm experiment 
such as the Square Kilometre Array (SKA) (with the pixel size chosen to be able to identify the wake signal, as 
described in the text) leads to a limit of G/i > 2.5 x 10~ 8 , an improvement over the current limit of Eq. p\ obtained 
from CMB anisotropy maps [3j. There are three important comments to make. First, a limit obtained by direct 
searches in position space is independent of some of the unknown parameters characterizing the scaling solution of 
strings (for example, it is largely independent of the number of strings per Hubble volume), whereas a bound obtained 
from the angular correlation function of CMB anisotropies depends on these parameters. Thus, a bound obtained 
by the methods we describe is more robust. Secondly with more statistics, even a signal strength below the noise 
can be detected. And finally, since cosmic string wakes produce signals with a distinguished geometry (namely thin 
wedges in redshift maps), these signals should be detectable by clever statistical techniques even if the amplitude is 
substantially smaller than the pixel noise. An example of a statistical tool which could be used to look for the specific 
string wake signal is the set of Minkowski functionals, statistics designed to characterize the topology of the structure 
of maps. For a preliminary study of this approach with references to the original literature see |24j . 



APPENDIX 



In this appendix we present how the (sin0)~ 2 factor in the expression for the wake brightness temperature arises 
from the line profile term 4>(is). We also explain why it does not imply any physical singularity in the limit 9 — > 0. 
We begin with the formula for the relative brightness temperature of photons (see pQ and references therein) 
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where the optical depth r„ is given by 
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(21) 



Here Njji is the column density of HI, and 4>{v) is the line profile. 

The line profile accounts for the nonzero width of the 21 cm line. The line profile is a sharply peaked function 
about ^21 =1420 MHz normalized such that 



4>{v)dv = 1 



(22) 



For the case at hand the line profile is not a Dirac delta function because the Hubble expansion leads to a velocity 
gradients along the line of sight. Due to this radial velocity gradient, the emitted 21 cm photons incur a local line 
broadening which one can interpret as a relative Doppler shift. 
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FIG. 5: A 21 cm light ray traverses a cosmic string wake of width w. 



In particular, let us consider photons reaching us at an angle 9 relative to the vertical of the wake having width 
w (see fig. Is!) . There is a relative Doppler shift 8v between photons from the front and back of the wake along the 
photon trajectory. The horizontal distance d along the wake between the back point of the photon trajectory and the 
photon trajectory exiting the wake is given by 



w tan 9 . 



(23) 



The relative Doppler shift is due to the Hubble expansion of the planar dimensions of the wake. Thus, the component 
of this expansion velocity in direction of the photon trajectory is given by 



v N 



Hd sin 



Hw sin 
cos 9 



(24) 



and we have a frequency width given by 



The line profile is 



8v = Vjj v<i\ (25) 



1 . 8v 8v . . 

— for vt i/ 2 i - ~7T ,V2i + ir , (26) 

0i/2 2 

and 0(i/) = otherwise. For a freq uen cy within the range of 21cm radiation passing through the wake the relative 



brightness temperature is given by (13). This makes the origin of the (sin#)~ 2 factor manifest. 

Any actual measurement will have a finite frequency resolution. If we integrate over a frequency interval larger 
than 8f, the (sin6>)~ 2 factor will cancel as it must. If we have a very small frequency resolution R, then the brightness 
temperature density at a fixed v will increase as 9 decreases until 8v becomes smaller than R. From that point 
onwards, the result will be independent of 9. Thus there is never any divergence associated with 9 — > 0. 
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